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Abstract 

Effect of two-particle correlations between impurities on the tem- 
perature of transition of a Fermi liquid in the superfluid or supercon- 
ducting state is analyzed. It is shown, that correlations with a radius, 
exceeding correlation length of the superconductor can have a pro- 
nounced effect on the transition temperature. The equation for the 
transition temperature is corrected for correlations of the impurities. 
Possible applications of the new equation are discussed. 

1. Effect of impurities or vacancies on thermodynamic properties of tradi- 
tional superconductors is well described by the Abrikosov and Gorkov (AG 
in what follows) theory of superconducting alloys [1J. Possible correlations 
among the impurities are neglected in this theory. Difference with respect to 
the pure superconductor in this theory is characterized by one parameter - 
electron mean free path I relative to the coherence length of the supercon- 
ductor £ • It has been observed recently [2] (and references therein), that 
certain impurities in the high-T c compounds demonstrate a tendency to 
short-range ordering, which influences the superconducting transition tem- 
perature T c of the compound. Even more clear example of such influence is 
presented by the superfluid 3 He in a high-porosity aerogel [31 0j. Aerogel is 
a rigid structure consisting of strands of a diameter a ~ 3 nm and the aver- 
age distance between the strands £ a ^40 nm for 98% porosity. The latter 
distance is usually taken as a characteristics of correlations within aerogel. 
Thermodynamic properties of 3 He in aerogel, in particular dependence of 
the suppression of the transition temperature with respect to that of pure 
3 He on £ /l demonstrate significant deviation from the universal behavior, 
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predicted by AG-theory. This deviation was attributed to the existence of 
structural correlations in aerogel with a characteristic radius which is compa- 
rable with £ [5] . In 98% aerogel depending on a pressure £o for 3 He varies 
from infinity at P m 5 bar to £o ~ 20 nm at solidification. The condition 
£a ~ £o is met at P m 15 bar. 

In the present paper effect of correlations on T c is considered under the 
assumption that the correlations are weak. Parameter, which characterizes 
deviation from the AG-theory turns out to be the ratio , where R is a 
length, characterizing correlation of impurities. For aerogel R ~ £ a and the 
effect of correlations on the T c can be significant. Analogous approach was 
used before for the situation when suppression of T c is small in comparison 
with the original T c [5J. Here this restriction is lifted and the argument is 
extended for arbitrary suppression of T c . This generalization is necessary 
for 3 He in aerogel since at a pressure about 6 bar the observed T c is zero 
so that the suppression of T c is equal to T c itself. 

2. Transition temperature T c is defined as a temperature of the onset 
of the long-range order. At this temperature appears a pole in the thermo- 
dynamic vertex part T(Q n , q) for scattering of two quasi-particles with the 
total momentum q = and the total frequency Q n . = [?]. Equation for 
T c has the standard form for all types of Cooper pairing: 

i = 11(0,0), (i) 

where A is a constant of the pairing interaction and polarization operator is 
defined as 

— -^(k+q, co m + n n )G(-k, u„ m ) =T^Lt ] (q, tK) , 

(2) 

G(k,u n ) etc. are thermodynamic Green functions of quasi-particles. In 
what follows L^(q, Q n ) always enters at Vl n = , and a shorthand notation 

L#?(q,0) = *m(q) is used. 

Effect of impurities is described by the potential 

(r) = 5>(r - r a ) (3) 

a 

which is a sum of individual potentials of impurities, situated in positions r a . 
These positions are random. In the standard technique [7] Green functions of 
quasiparticles are expanded in the Born series over U(r) and averaged over 
positions of the impurities. The sum of the obtained series can be written 
as: 

(C?(e;k,kO) = 27r5(k-k')[G 1 - (E)]" 1 , (4) 
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Figure 1: 



where (£) is the averaged self-energy. Two first terms of the Born series 
for (E) are represented by the diagrams Fig. 1. Correlations among the 
impurities introduce changes in the procedure of averaging. 

To every cross corresponds a function u(k!— k) Y, a e* < - k ~ k ') ra , where k and 
k' - the in-coming and out-going momenta. Essential contribution comes 
from the term with two crosses. It is proportional to the averaged sum: 
(E«,i) e! ' k ~ kl ' r ° +i ' kl ~ k ' r ')' Summation over r a renders a factor (2ir) 3 n5(k — 
k') . The remaining sum is a structure factor 



where q = kx — k. For uncorrelated impurities S(q) = 1 + (27r) 3 n<5(q) . 
The unity comes from the summand with r& = r a . The term with 5(q) can 
be dropped, because in further calculations it enters with the factor u(0) , 
which is excluded by a shift of the ground state energy. When correlations are 
present there appears an additional term 5(q) which is directly related to the 
correlation function in the coordinate representation. Summation in Eq.(5) 
for Tf, 7^ r a can be changed for integration with the probability w(rb\r a ) 
to find a particle in the point rj, if there is a particle in the point r a . In 
the isotropic case w(rb\r a ) depends on a relative distance r = |r& — r a | . At 
r — > oo correlations vanish and w(r) tends to a constant. Normalization of 
w (r) is chosen so that this constant is unity. A measure of correlations is 
v(r) = w(r) — 1 (cf. [8J). With these notations 



S 2 (e; k) = n / u(k a - k)u(k - k x )[l + S(k - ki)]G?(e; k 1 )— ^. (7) 



Integration of two crosses in one block £2 is graphically represented by the 
dashed line, connecting the crosses (cf. Fig.l). To this line corresponds now 



S(q) = (^ e iq(r "- rtl) ) 



(5) 



b 




(6) 



and: 
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the sum 1 + S(q) . The averaged Green function has the same form as for 
uncorrelated impurities 

(G(w n ;k,k')) = (2n) 3 S(k-k , )(iuj n -pi + ^sgnw n )-\ (8) 

It 

except that the average inverse time between the collisions includes S(q) : 

- = 2nnu I \ U (9)\ 2 [1 + S(9)]^. (9) 

T J 47T 

Here 9 is the angle between ki and k .The new term S{9) takes into ac- 
count effect of interference of de Broigle waves of quasiparticles scattered by 
different impurities. This effect can be formally included in the definition of 
cross-section |w(#)| 2 — l M (#)| 2 [l + S{9)\ . It is determined by the properties 
of the normal phase and does not depend on T c . The effect of interference 
is usually small, since the average distance between the impurities a is large 
in comparison with the wave length of quasiparticles, and 5(q) vanishes at 
|q| > 1/a. Contribution of S'(q) to the integrand in Eq. (7) is comparable 
with unity only at angles 9 ~ 1/akp . For small concentration of impurities 
x the relative contribution of correlations is on the order of x 2//3 . 

At the averaging of the polarization operator except for substitution of 
the averaged Green functions in Eq.(2) it is necessary to take into account the 
contribution of "ladder" diagrams of the type represented by the first diagram 
in the r.h.s. of the equation Fig.3. The dashed lines bring in the factor 
1 + S(k — k^ , which can also be included in the cross-section |w(ki — k)| 2 . 
For definiteness the explicit analysis here is made for the p-wave pairing. 
That makes possible to apply the obtained results to the superfluid 3 He in 
aerogel. The interaction leading to the p-wave pairing is proportional to the 
scalar product of unit vectors k, k/ : V(k, k/) = 3Vi(k-k/) . As the result the 
loops on the left and the right ends of the ladder diagrams acquire factors 
k and k/ respectively. The cross-section of scattering on the spherically 
symmetric impurities can be expanded in a series of Legendre polynomials: 

Kkx-k)! 2 = a + 3<7i(k • k/)) + 5a 2 P 2 (k ■ kf)) + .... (10) 

After the angular integration in every loop of the ladder diagrams only 
o\ -component remains. It enters results via inverse collision time 1/ri = 
27rz/o0"i • Summation of the series for n(0, 0) renders familiar equation for 
T c at the p-wave Cooper pairing: 

l = 27rn,ViV 1 , (11) 

n > W„ + i./ZT tr 
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Figure 2: 



where l/r tr = l-nnv^a^ — ax) . For the pairing with other angular momenta 
I 7^ 1 a proper o\ and I/77 have to be used. In particular for the case of 
s-wave pairing the resulting 1/r would be zero and transition temperature 
is not lowered in agreement with Refs. [TJ [U] 

So, within the scheme of the AG theory, correlations among the impurities 
enter the equation for T c only implicitly, via the mean free path, determined 
by the normal phase. 

3. A non-trivial effect of correlations occurs when the contribution of 
diagrams with four crosses is taken into account in Il(q, Q n ) . Of particular 
interest for the present discussion are the diagrams with three dashed lines 
as shown on the Fig. 2. 

Without correlations it would correspond to the fourth order correction to 
the amplitude of scattering of two quasiparticles with zero total momentum 
on the same impurity. For interacting quasiparticles in the intermediate 
state instead of the product of two single particle Green functions G n (ki + 
qi)GL n (— ki) the two-particle Green function G H (ki + qi, — ki,co>;, —cji] k 2 + 
qi, — k2, co m , — co m ) has to be used. As a function of qi G u is singular at 
qi = . This singularity corresponds to formation of a virtual Cooper pair. 
It increases the weight of states with q! ~ l/£ in the integral, corresponding 
to this diagram. If the impurities are correlated on a distance ~ £ there 
remains a net effect of interference of the waves with qi ~ l/£o scattered by 
impurities. Depending on a sign of correlations it can decrease or increase 
the destructive effect of impurities. The aim of the following discussion is to 
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Figure 3: 

find out when the effect of interference is essential. 

Interaction of electrons conserves the total momentum q! and the total 
frequency Q = but, unlike the elastic scattering by impurities, it does not 
conserve each of the frequencies ui . As the result the polarization operator 
became a sum over two indices: 11(0) = T 2 J2im Li m (0) , where L; m (0) is a 
block, represented by the black rectangle in the Fig. 3 . 

In analytic form the equation Fig. 3 can be written as: 

L m = U0) -t^T^MmjLj, (12) 

j 

where 

M m3 = no*J Uq)r(q)[l + ^(q^q)^, (13) 

Lj = TJ2 m L 3m (0) and t,(q) = Bj(q)/[1 - rmB^q)] with 5,(q) = 
^ arct 9 '(2^1) • In a limit 9 ^ it tends to tj(0) = 27ri/ /|a>jj with 
Cjj = Uj + l/2r tr . 

The matrix M m j is proportional to a\ , which brings in a small factor 
(Co/hr) 2 ■ A special situation arises when the impurities are correlated on 
a distance R which is greater than £ . In this case S(q) is enhanced 
in a region q ~ 1/R . In the limiting case R ^> £ an functions under 
the integral in Eq.(13) except for S(q) and T(q) can be taken at q — 
so that M mj = na 2 t m (0)Qt i (0) , where Q = jr(q)S(q)A . At the 
transition temperature and at small q T(q) = —7/^0 (Co?) 2 with coefficient 
7 = 12/7^(3) . A singular behavior of T(q) at q — > increases the weight 
of small q in the integral and since S(q) is enhanced at q ~ 1/R the value 
of the integral is also enhanced. For estimation of integrals, containing S 
the form S(q) = A5(q — go) with qo — 1/R is convenient. The coefficient 
A is determined from the normalization condition / S(q)d 3 q = (27r) 3 nv(0) , 
so that A = 2ir 2 R 2 nv(0) . Omitting insignificant coefficients on the order of 
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unity we arrive at 

fl 2 m,(0) 

(14) 

This estimation shows, that the second term in the rhs of Eq.(12) relative 
to the first is on the order of (R/l) 2 ■ At i? <C / Eq.(12) can be solved by 
iterations. As a first iteration we have: 

L m = t m (0) - t m (0)T]T M m ^(0), (15) 

3 

Summation of this equation over m with account of Eq.(l) after standard 
transformations renders the equation for for T c : 

T AG " V2 ' 47rT AG r tr 



where 

r 2 



f(T a ) = [ 

Vq J 



T£t m (0)t m (q) 



V(q)[l+S(q)]^. (17) 



and ip(z) = ^lnT(z) is di-gamma function. Equations (16) and (17) form 
a closed system, which for given material parameters and a given form of 
the structure factor determines the transition temperature in a liquid with 
correlated impurities. Difference with respect to the standard result Tag is 
represented by the function f(T c ) in the r.h.s. of Eq.(16). 

Further iterations of Eq. (12) render consecutive terms in the expansion 
of the solution of Eq. (12) in powers of the parameter (R/l) 2 ■ The accuracy 
of Eq. (12) itself is limited by a different parameter. The kernel (13) takes 
into account interaction of the virtual Cooper pair only with two correlated 
impurities. Averaging of diagrams including interaction with three correlated 
impurities requires knowledge of the three-particle correlation function. For 
estimation of the contribution of this process we can decouple the three parti- 
cle correlation function in a product of the two-particle correlation functions 
and keep only the most singular terms. The correction would contain extra 
factor Q£o/kr ~ R 2 /(,okr ■ Higher order correlations can be neglected if this 
parameter is small. 

For unconventional superconductors the actual situation is £ < hr ■ We 
assume that the strong inequality £ ^ hr is met. Then there exist a window 
£o R <^ kr ■ In a limit R ^> £o the asymptotic expression for f(T c ) at 
(R/Co) oo can be used. The function t m (q) in the integrand of Eq. (17) 
can be taken at q = and the sum over m can be done explicitly 

r ?^°'=I^G + sk) < 18 > 
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The remaining integral over q in the limit R ^> £o is proportional to 
(R/2li) 2 with a coefficient a ~ 1 depending on the explicit form of the 
structure factor. The resulting form of f(T c ) is then: 



The sign of a depends on a limiting value of the correlation function v(r) 
at r — > . At v (0) > a < , i.e. a tendency of impurities to form groups 
favors formation of Cooper pairs and increases T c and vice versa. 

4. For application of the obtained formulae to aerogel we have to sub- 
stitute its structure factor in the definition of Q . 98% Aerogel within the 
interval of distances a <C r <C R has a fractal structure with the dimen- 
sionality D — 1.6 -T- 1.8 [5]. Its structure factor within 1/R <C q <C 1/a 
varies as q~ D . As a consequence the integral Q ~ / dq/q D is dominated by 
the lower limit ~ 1/R. Unfortunately, the available data for the structure 
factor are given in arbitrary units so that the value of Q can be found only 
up to unknown factor. Instead of the measured structure factor a model 
expression was suggested j6], assuming that at r > R correlations decay 
as exp(— r/R) , R in this model is a fitting parameter. Simplified analysis, 
using the limiting expression valid for R ^> £ can be used only at pressures 
above 20 bar. In this region Eqns. (16), (19) reproduce qualitative features 
of the experimentally observed dependence of T c on pressure. In particular, 
the difference T c — Tag is practically constant within this interval, but the 
interval itself is not large enough for making definitive statements. In a re- 
gion P < 20 bar the limiting expression for t m (0) is not sufficiently accurate 
and Eqns. (16), (17) have to be solved with the full t m (q) . That requires 
rather involved numerical calculations which are not done yet. 

In conclusion, it is shown that effect of correlations on the suppression of 
the transition temperature of a Fermi liquid in the superfluid (superconduct- 
ing) state by impurities can be essential if correlation radius R of impurities 
is greater than the correlation length of a superconductor £0 • Correlations 
can increase the transition temperature with respect to the value, given by 
the standard theory of superconducting alloys. The obtained equations ex- 
press the transition temperature in terms of the structure factor of ensemble 
of impurities. The quantitative analysis takes into account only two-particle 
correlations. Within this approximation a relative change of T c is character- 
ized by the parameter R 2 /^ol tr . This parameter simultaneously determines 
accuracy of the approximation. If R 2 /C,ol tr ~ 1 , the results can be used only 
for qualitative estimations. Good candidates for observation of the discussed 
effect are superconductors with small values of £ • 




(19) 
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